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This  paper  is  concerned  with  eigenvalue  problems  for  boundary  value 


problems  of  ordinary  differential  equations  posed  on  an  infinite  interval. 
Problems  of  that  kind  occur  for  example  in  fluid  mechanics  when  the  stability 
of  laminar  flows  is  investigated.  Characterizations  of  eigenvalues  and 


spectral  subspaces  are  given  and  the  convergence  of  approximating  problems 
which  are  derived  by  reducing  the  infinite  interval  to  a  finite  but  large  one 
and  by  imposing  additional  boundary  conditions  at  the  far  end  is  proved. 


Exponential  convergence  is  shown  for  a  large  class  of  problems. 


AMS(MOS)  Subject  Classification:  34B05,  34B25,  34C11,  34E05,  34A45 

Key  words:  Boundary  value  problems  of  linear  equations,  spectral  theory 
of  boundary  value  problems,  boundedness  of  solutions, 
asymptotic  expansion,  theoretical  approximation  of  solutions 

Work  Unit  No.  3  -  Numerical  Analysis  and  Computer  Sciences. 


Sponsored  by  the  United  States  Army  under  Contract  No.  DAAG29-80-C-0041 .  This 
material  is  based  upon  work  supported  by  the  National  Science  Foundation  under 
Grant  No.  MCS-7927062. 


SIGNIFICANCE  AND  EXPLANATION 


This  paper  deals  with  eigenvalue  problems  for  ordinary  differential 
equations  posed  on  an  infinite  interval.  These  problems  have  the  following 
form:  We  have  a  system  of  linear  ordinary  differential  equations  depending 
linearly  on  an  (eigen)  parameter  and  boundary  conditions  at  a  finite  point 
tQ,  and  we  look  for  nontrivial  solutions  which  fulfill  the  differential 
equation  on  the  interval  [tQ,“)  and  the  boundary  condition  at  tQ,  and 
which  have  a  finite  limit  at  infinity.  Such  problems  occur  frequently  in 
fluid  mechanics  when  the  stability  of  flows  over  infinite  media  is 
investigated.  In  this  paper  eigenvalues  and  spectral  subspaces  are 
investigated  and  characterized.  A  suitable  way  to  solve  such  problems 
numerically  seems  to  be  to  cut  the  infinite  interval  at  a  finite,  large  enough 
point  T  and  to  impose  suitable  'artifical'  boundary  conditions  at  t  =  T  in 
order  to  get  an  eigenvalue  problem  for  a  two-point  boundary  value  problem  on 
[t0,T] .  The  questions  that  arise  immediately  are:  What  boundary  condition 
at  T  achieves  convergence  of  the  'finite'  solution  to  the  'infinite' 
solution  as  T  converges  to  infinity,  and  what  order  of  convergence  can  we 
expect.  In  this  paper  exponential  convergence  is  shown  for  well  constructed 
approximating  problems. 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
summary  lies  with  MRC,  and  not  with  the  author  of  this  report. 


(1.3)  y  £  C((1,«)  )  :  <“>  y  £  C([1,"))  and  lim  y(t)  exists 

where  the  nxn  matrices  A,G  £  C( [1 »“] )  and  A(“)  *  0  .  A  sound  theory  for  inhomogenous 
boundary  value  problems  on  infinite  intervals  has  been  developed  (3ee  Lentini  and  Keller 
(1980),  de  Hoog  and  Weiss  ( 1980a, b),  Markowich  ( 1980a, b,c))  but  not  much  attention  has  been 
paid  to  eigenvalue  problems  with  a  singularity  of  the  second  kind,  de  Hoog  and  Weiss 
(1976)  established  a  theory  for  eigenvalue  problems  in  the  case  that  the  differential 
equation  has  a  singularity  of  the  first  kind  (a  *  -1)  and  that  G(“>)  *  0.  They  could  show 
that  the  spectrum  has  no  finite  limit  point  and  that  the  spectral  subspaces  associated  with 
a  particular  eigenvalue  are  finite  dimensional.  They  also  considered  difference  schemes 
for  problems  which  have  been  transformed  to  a  finite  interval,  and  they  derived  convergence 
results  for  eigenvalues  and  spectral  subspaces  U3ing  the  collective  compactness  of  the 
difference  schemes.  They  also  derived  certain  properties  of  the  spectrum  and  the  spectral 
subspaces  of  (1.1),  (1.2),  (1.3)  under  the  assumption  that  the  matrix  A(“>)  +  1G(<“)  has  no 
eigenvalue  on  the  imaginary  axis  (see  de  Hoog  and  Weiss  (1980a)).  They  showed  that  all 
eigenvalues  \  for  which  this  assumption  on  A(“>)  +  1G(“)  holds  are  isolated  and  their 
spectral  subspaces  are  finite  dimensional.  Their  proofs  hinge  on  the  Fredholm  property  of 
the  differential  operator. 

Sponsored  by  the  United  States  Army  under  Contract  No.  DAAG29-80-C-004 1 .  This 

material  is  based  upon  work  supported  by  the  National  Science  Foundation  under 
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The  goal  of  this  paper  Is  twofold.  First  to  derive  properties  of  the  spectrum  and  the 


generalized  eigenvectors  of  (1.1),  (1.2)/  (1.3)  and  second  to  consider  the  approximating 
eigenvalue  problems 

(1.4)  -  taA(t))^,  -  ,  1  <  t  <  T,  T  >  1 

(1.5)  BXjd)  -  0 

(1.6)  smx^T)  -  0  . 

These  problems,  with  a  suitably  chosen  matrix  S(T)  are  'regular*  two  point  boundary 
eigenvalue  problems  which  can  be  solved  by  any  appropriate  code,  for  example  by  collocation 
(see  de  Boor  and  Swartz  (1980)).  A  class  of  matrices  S(T)  for  which  the  eigenvalues  and 
spectral  subspaces  of  (1.4),  (1.5),  (1.6)  converge  to  those  of  (1.1),  (1.2),  (1.3)  will  be 
defined  and  the  order  of  convergence,  which  turns  out  to  be  exponential  in  the  most 
important  cases,  will  be  estimated. 

This  paper  is  organized  as  follows.  In  Chapter  two  the  case  where  A(°°)  has  no 
eigenvalue  on  the  imaginary  axis  and  where  G(“)  -  0  is  treated,  in  Chapter  three  the 
assumption  G(“)  «  0  is  neglected.  In  Chapter  four  no  assumption  on  the  eigenvalues  of 
A(")  are  made,  but  a  certain  order  of  convergence  of  G(t)  to  0  is  required.  In 
Chapter  5  the  Orr-Sommerfeld  equation,  a  fluid  dynamical  problem  posed  as  an  eigenvalue 
problem  on  an  infinite  interval,  is  dealt  with  and  appropriate  approximating  problems  are 
devised. 

It  is  of  particular  interest  that  the  approximation  theory  in  the  case  G(*“)  *  0  is 
treated  by  using  Griegorief f ' s  (1975)  'discrete'  approximation  theory  for  eigenvalue 
problems,  which  allows  the  approximating  operators  to  be  defined  on  different  spaces  which 
-  in  some  sense  -  converge  to  the  space  on  which  the  eigenvalue  problem  is  posed.  This  way 
of  pursuing  simplfies  the  analysis  essentially. 


I 


2.  G(»)  -  0  i  The  'Compact'  Case 

We  assume  that  A(»)  11m  A(t) 

t»» 

G(“)  :■  lim  G( t)  “  0  .  We  transform 
t+® 

(2.1) 


has  no  eigenvalue  with  real  part  zero  and  that 
At*)  to  its  Jordan  canonical  form  j(“) 

A(-)  -  FJ(-)F_1 


and  assume  that  J («•)  has  the  blocks  structure 

(2.2)  J(«)  -  dlag(J*  ,  jj 

where  J*  contains  all  Jordan  blocks  which  have  eigenvalues  with  real  part  larger  than 
zero  and  contain  all  Jordan  blocks  with  eigenvalues  with  negative  real  part.  Let 

J*  be  a  r+  *  r+  matrix  and  J^  a  r_  *  r_  matrix  and  let  D+  resp  D_  be  the 
projection  onto  the  sum  of  invariant  subepaces  associated  with  the  eigenvalues  of 
J*  resp  J^  .  We  define  a  solution  operator  h  of  the  problem 

(2.3)  z'  -  t°J(»)z  +  t“g(t) 


for  all  real  a  >  -1  as  follows 


(2.4)  (Hg)  ( t) 

where  6  >  1  and 


♦  (t)  D  ♦  1  (  s)  a"’g(  s)ds  +  $<t)  ft  D  ♦  1  ( s)  sQg(  e)  ds 

-  i 


(2.5)  *<t)  -  «p(^  t“+1)  . 

This  operator  has  been  used  by  de  Hoog  and  Weiss  ( 1980a, b)  and  they  showed  that 
Hg  e  C (  ( 1  ,<“]  )  if  g  e  C( [ 1 ,®] )  and  that 

(2.6)  <Hg!<«)  -  -J(»)-1g(«) 

holds.  Moreover  H  s  C(  (£,*])  ♦  C((6,®1)  is  bounded  and  its  norm  is  Independent  of  6 
and 


(2.7) 


I  (Hg)  (t)  I  <  const. (Igl  t  +  lexp(^  f2^'1 )  t0*-*" 1  j  I  Igl  ) 

[j,®]  2 


holds  for  t  >  2.  As  norm  in  C(  (£,"])  we  take  1*1  .  ,  which  denotes  the  max-norm  on 

l  o  ,®  J 


the  interval  [S,®]  . 

Now  we  investigate  the  problem 


(2.8) 

(2.9) 


y'  -  t°A(t)y  -  taG(t)f(t)  1  <  t  <  - 


By< 1 )  «  0 


-3- 


(2.10) 


y  e  C< [1,»] ) 


where  f  e  C([1,“l).  Substituting  u  -  F_1y  the  general  solution  of  the  transformed 
problem  (2. 8),  (2. 10)  is 


(2. 1 1 ) ( a) 


u ( t )  -  $(t)G  5  +  (H(J  -  J(“))u)(t)  +  (HF  Gf)(t),  See 


where  J(t)  =  F_,A(t)F  and  the  n  *  r_  -matrix  G_  is  obtained  from  D_  by  cancelling 
all  columns  which  have  only  zero  entries. 

Obviously  the  operator 


( 2 . 1 1 )  ( b) 


H(J  -  J(“))  :  C<  («,“])  +  C([6,“)) 


fulfills  IH(J  -  J(«))ljj  <  1  for  $  sufficiently  large.  Therefore  u  in  (2.11)  is 

defined  uniquely  on  [S,“]  and  can  be  extended  uniquely  to  11,“]  . 


Setting 


(2.12)  <(i_(t)  =  ((I  -  H (J  -  J(“))  $G_)(t) 

(2.13)  *(Gf)(t)  -  ((I  -  H(J  -  a  {<»>))  ~ 1  HF~ 1  Gf )  ( t ) 
we  get  the  general  solution  of  (2.8),  (2.9),  (2.10): 

r 

(2.14)  y(t>  =  Fi|/_(t)C  +  Ft|i(Gf)(t),  t  £  [1,“1,  icC~. 

So  (2.8),  (2.9),  (2.10)  is  uniquely  soluble  for  every  f  e  C([1,“I)  if  and  only  if  the 


r  x  r  matrix 


BFi)i  (1)  is  regular. 


B  is  assumed  to  be  a  r  x  n  matrix.  (2.6)  and  (2.11)(a)  imply  that  y(“)  =  0  . 
Now  we  define  the  operator  V  as  follows: 


C(  (1 ,“]  )  -►  C(  [1,“]  ) 


f  ♦  Vf  «  y 


where  y  is  the  solution  of  (2.8),  (2.9),  (2.10).  V  is  defined  properly  if  and  only  if 

(2.15)  holds.  This  is  no  restriction  because  if  X  =  0  is  an  eigenvalue  of  (1.1),  (1.2), 
(1.3)  we  substitute  X  »  X  +  y  so  that  the  problem  with  Aft)  replaced  by  Alt)  +  XG(t) 
has  not  y  -  0  as  eigenvalue. 

(2.7)  and  (2.14)  imply  that  V  is  bounded. 


r- 


Obviously  the  eigenvalue  problem  (1.1),  (1.2),  (1.3)  is  equivalent  to 


(2.17)  Vf  -  uf 
with 

(2.18)  w  “  X  • 

Our  goal  is  to  show  that  V  is  compact.  We  need  the  following 

Lemma  2.1:  Assume  oeC([S,~]).o(t)-*0  as  t  ♦  o  >  0  and  S  >  1.  then  the  set 
Aq  defined  by 

A  i-  {f  e  C(  (1,-1  )  |lfl  .  <  C,,  If(t)l  <  C,o(t>  for  t  >  6,  If  ( t)  I  <  C.t“}, 
a  [1,"]  1  2  3 

for  aeR  is  conditionally  compact  in  C([1,“])  . 

Proof.  Given  e  >  0  we  choose  T  “  T(e)  >  5  so  large  that  oft)  <  —  for  t  >  T  . 

2 

Obviously  there  is  a  finite  collection  of  intervals  1^  for  i  *  1(1)(N(c)  -  1)  ,  whose 
conjunction  is  [1,T],  and  there  are  points  in  lj  so  that 

sup  sup  lf(t.)  -  f(t)l  <  e  ,  i  «  1(1  >(N(e)  -  1)  . 

fcA  tel  ,  1 

o  i 


This  is  fulfilled  if  1 |  <  E/fCjT01)  with  tj^  arbitrary  in  . 

Setting  t  ■  “  Theorem  5  in  Dunford  and  Schwarts  (1957)  is  fulfilled  and  the 
N(  e ) 

Lemma  follows. 

From  (2.11)  and  (2.7)  we  conclude  that 


I (Vf ) ( t) I  <  const ( I $( t)G  I  +  IJ(t)  -  J(-)l  +  IG( t) * 

if  if 


t  >  6  >  2 


holds  because 


in  (2 •  1 1 ) ( a)  equals  -(BE\|»  (1))  E\(/(f)(1) 


Setting 


(2.20)  oft)  -  l$(t)G  I  +  IJ(t)  -  J(-)l  +  lG(t)l  +  »  exp(-^-(-  1  )  t3* 1 )  I 

if,-)  if,-]  2*+1 

we  notice  that  o(t)  +  0  as  t  ♦  -  and  therefore 

(2.21)  { Vf | f  e  C([1,-]>,  lfl(1  <  1}  £  Aq 

for  some  constants  C2,  C3.  So  V  is  a  compact  operator  on  C([1,«]). 


r 


wm 


Because  0(11,°°])  is  infinite  dimensional  and  V  is  bounded  the  spectrum  o(V) 
consists  of  an  infinite  sequence  of  eigenvalues  <  0  of  finite  algebraic  multiplicities 
and  v  •  0  e  o(V)  ia  the  only  accumulation  point  of  the  u^'s  •  This  implies  that  the 
eigenvalues  1  *  —  of  (1.1),  (1.2),  (1.3)  have  no  finite  accumulation  point,  they 


fulfill 

(2.22)  11^1  ♦  •>  as  i  ■»  «•  . 

The  spectrum  of  compact  operators  is  described  in  Dunford  and  Schwartz  (1957),  Chapter  VII, 
Theorem  5. 

Let  u  *  0  be  a  fixed  eigenvalue  of  V  .  We  want  to  investigate  the  spectral 
subspace  associated  with  u  .  The  spectral  projection  is  given  by 

(2.23)  E  -  E(  U )  -  ——  I  (z  -  V)-1dz:C(  [1,»]  ♦  C[{1,-)> 

2n  X  j. 


where  f  is  a  circle  centered  at  u  which  contains  no  other  eigenvalue  of  V  .  Moreover 

(2.24)  rank(E(U) )  -  m 
where  m  is  the  algebraic  multiplicity  of  u  .  Let 

g 

(2.25)  Range  (E)  »  spant  *,*••,  *P  }  “  N((u  -  V)  ) 

i  cn 

hold,  where  the  are  generalized  eigenfunctions  of  V  .  N  denotes  the  null  space  and 

8  the  ascent  of  u  -  V  . 

As  the  range  of  E  is  invariant  under  V  we  get 

m 

(2.26)  *>'  -  taA(t)V>  =»  l  a  taG(t)*  ,  Be  ( 1  )  -  0  ,  <  E  C(  [  1  . 

1  1  j»1  l  1  1 

The  m  x  m  matrix  (a^j)  can  156  assumed  to  be  in  Jordan  canonical  form  with  the  only 

eigenvalue  X  =  ^  .  This  can  always  be  achieved  by  a  basis  transformation.  So  every 

element  'f.  is  contained  in  a  finite  chain  ^  .•"»  •f  which  fulfills 

ri  ri 

(2.27)  -  t“(A(t)  +  XG(t>)*  -=  0  ,  Bp  (1)  «  0  ,  *  E  C((1,»)) 

ri  ri  ri  ri 


(2.28) 


(2.29) 


i f *  -  t“(A(t)  +  XG<  t )  )gj  -  taG(t)^  ,  (1)  =  0  ,  <f  E  C(11,“)) 

r2  r2  ri  r2  r2 

•  •  * 


<?•  -  ta(A(t)  +  XG(  t)  )*fi  -  t°G(t)^  ,  &p  (1)  >  0,  *  t  C(  (1,«]  )  . 

ri  ri  *1-1  Tl 


I 


t 


-6- 


Obviously  every  ^  fulfills 

(2.30)  I  #  (t)l  <  C^ft) 

but  a  stronger  estimate  will  be  advantageous  for  the  order  of  convergence  of  the 
approximating  eigenvalue  problems  (1.4),  (1.5),  (1.6).  Therefore  we  assume  that 

( a+ 1 ) k  +1 

(2.31)  a(  7  ),  g(  7  )  e  c  ([0,  ■£])  ,  &  >  1  ,  a  e  NQ 

holds,  where  k_  is  the  largest  algebraic  multiplicity  of  an  eigenvalue  of  A(“)  with 
negative  real  part.  Under  this  assumption  it  follows  from  Markowich  (1980b)  that  4 

1 

which  is  the  solution  of  the  homogenous  problem  (2.27)  as  well  as  the  4  for  i  >  1 

i 

which  are  solutions  of  inhomogenous  problems  decay  exponentially  and  they  fulfill  the 
estimate 

(v  +  e) 

(2.32)  l*r  ( t ) I  <  C1  exp( - -  t“+1),  t  >  6 

where  v_  is  the  largest  negative  real  part  of  eigenvalues  of  A(“)  and 

t  -  e(5)  >  0  fulfills  v  +  e  <  0  and  e  ( 6 )  •»  0  as  S  ♦  -  .  Therefore  all  elements  in 
Range  (E(t))  fulfill  the  estimate  (2.32). 

Now  we  want  to  investigate  the  convergence  of  the  eigenvalue  and  generalized 
eigenvectors  of  the  approximating  problems  (1.4),  (1.5),  (1.6).  As  a  notion  of  the 
distance  of  closed  subspaces  we  use  the  'gap'  (see  Osborn  (1975))  which  is  defined  as 
follows 

(2.33)  gap(M,N)  «  max(  sup  dist(x,N),  sup  dist(M,y)) 

x  £  M  y  c  N 

1x1-1  lyl-1 

where  M,N  are  closed  subspaces  of  a  Banach  space  (X,*,l)  and  diet  is  defined  as 

(2.34)  dist(x,N)  -  inf  Ix-yl  . 

y£N 


We  define  the  operators  V_  for  T  sufficiently  large  by 


(2.35) 


V„  1 
T 


C< (1 ,“] )  *  C( 11,“)) 


* 

T  T 


where  xT  fulfills 


(2.36) 


X^,  -  t'VtlXy  -  taG(t)f(t),  1  <  t  <  T 

(2.37)  BXjd)  -  0 

(2.38)  S(TJXj,(T)  -  0 
and 


(2.39)  x^tj  -  x^T)  for  t  >  T  . 

This  definition  makes  sense  if  and  only  if  (2.36),  (2.37),  (2.38)  is  soluble  for  every 
f  £  C( [1,-J )  and  T  sufficiently  large,  de  Hoog  and  Weiss  (1980b)  have  shown  that  this 
is  the  case  if  (2.15)  holds  and  the  r  «  n  -matrix  S(T)  fulfills 

(2.40)  IS(T)K  <  const,  as  T  *  “ 

(2.41)  l(S(T)FG+)  %  <  const  as  T  *  * 

where  the  n  x  r+  matrix  G+  is  obtained  by  cancelling  all  columns  of  D+  which  have 
only  zero  entries.  Moreover  the  stability  estimate 

(2.42)  *V[1,T]  *  COBSt*<,f' [1<T]  +  'Y(T),> 

holds  for  problems  of  the  form  (2.36),  (2.37)  and 

(2.43)  S(T)xl<T)  =  y(T> 

instead  of  the  homogenous  boundary  condition  (2.39).  de  Hoog  and  Weiss  (1980b)  have  also 
shown  that  (2.41)  is  necessary  if  (2.40)  holds  and  they  constructed  matrices  S(T) 
fulfilling  (2.40),  (2.41)  more  explicitly.  Obviously  the  estimate  (2.42)  with  y(T)  =  0 
and  the  definition  of  VT  imply 


(2.44) 

Every  operator 

(2.45) 

holds  and  V^f 
(2.8),  (2.9), 

(2.46) 

(2.47) 


•  V  I  r .  ,  <  const . 

T  11,"] 

vT  is  compact  because 

1 ( V  f ) 1 1  <  const.  T°lf«r, 

T  [1 ,T)  [ 1 ,T] 

is  constant  on  [T,“)  .  By  adding  the  identity  S(T)y(T)  =  S(T)y(T)  to 
(2.10)  and  by  subtracting  from  (2.36),  (2.37),  (2.38)  we  get  the  problem 


<*T  "  y)  '  "  t  A(t)(xT  -  y)  «  0  ,  1  <  t  <  T 

B(XT  -  y)  ( 1 )  -  0 
S(T)(xt  -  y) (T)  =  -SlT)y(T)  . 


(2.48) 


Applying  estimate  (2.42)  implies 


IV  f  -  Vfl,,  ,  <  const. ly(T)l  <  const  o(T)lfl„ 

T  [1,T]  (1,° 


(2.49) 

where  c  is  defined  in  (2.20).  Also  we  get 

(2.50) 


IV  f  -  Vfl ,, 

T  [1,“ 


<  21  Vf  -  Vfl  ,,  +  2 1 y I  .  , 

T  [1,T]  J  (T,“l 


because  of  (2.39).  Therefore  V,y  converges  to  V  in  the  norm  and 


VI , .  ,  <  const .C(T) 

[  1  ,°°1 


V  f  *  uf 
T  T 


1 


(2.51)  IV 

T 

holds. 

It  should  be  noticed  that  G(t)  *0  as  t  +  “  is  absolutely  crucial  for  the  norm 
convergence. 

The  eigenvalue  problem  (1.4),  (1.5),  (1.6)  is  equivalent  to 

(2.52) 
with 

(2.53)  u 

i  ^  _ 

T 

The  generalized  eigenfunction  of  (1.4),  ('.5!,  (1.6)  are  obtained  by  restricting  the 
generalized  eigenfunction  of  (2.52)  to  (1,T). 

Because  of  the  compactness  of  VT  there  is  an  infinite  sequence  of  eigenvalues 
*  0  accumulating  at  0  .  The  compactness  and  norm  convergence  (2.51)  allows  us  to 
apply  Osborn's  (1975)  result.  We  derive  the  following  convergence  statements. 

For  T  sufficiently  large  there  are  exactly  m  eigenvalues  of  VT  ,  counted 
according  to  algebraic  multiplicities,  which  lie  in  the  circle  T  defined  in  (2.23). 
These  m  eigenvalues  converge  to  u  as  T  *  "  .  The  spectral  projections 


(2.54) 
fulfill 

(2.55) 

Setting 

(2.56) 
we  get 

(2.57) 

where  (V„  -  V) 


h  /  (z  ■  VJ  1<3z 


2ni 


(a)  rank(ET) 


rank(E) 


(b)  gap(ET,E) 


0  as  T  ♦ 


-  i  < 

»  v  i 

li  *  —  ;  u 

T  m  >  *T 

1=1 


max( luT  -  ul,  gap(ET,E))  <  const. I (VT 
Range(E)  denotes  the  restriction  of  VT 


V) 


Range(E)  (I,®) 
V  to  Range(E). 


Resubstituting  X  -  —  ,  and  assuming  that  (2.31)  holds  we  get  the  following 


estimate  which  implies  exponential  convergence: 


(v  *e) 


(2.58) 


max( 


|X  -  ^r—  |,  gap  (E^.E) )  <  const.  exp(-^ —  T°+1)  . 


The  estimate  for  the  order  of  convergence  of  each  X*  to  X  is  worse: 


(2.59) 


i  (V  +  £) 

IX  -  XTI  <  const.  exp(-;-+lTr 


Q+  1  \ 

T  J  ,  i  »  1 ( 1 )m  . 


Here  0  is  again  the  ascent  of  u  -  V  .  The  constants  in  (2.58),  (2.59)  are 
independent  of  T  but  may  very  well  depend  on  X  .  If  the  assumption  (2.31)  is  dropped 
it  is  easily  derived  from  Markowich  (1980b),  Chapter  2  that  the  generalized  eigenfunction 
decay  faster  than  every  (negative)  power  of  t  and  so  the  same  is  true  for  the  order  of 
convergence  (2.58),  (2.59). 

A  possible  choice  for  s(T)  is 
(2.60)  S(T)  5  S  =  (G+)TF_1 

where  the  superscript  T  denotes  transposition-  The  condition  (2.41)  is  fulfilled  because 
(2.61  )  SFG+  *  Ir 

+ 

holds  for  the  choice  (2.60),  which  has  been  used  by  de  Hoog  and  Weiss  (1980b)  for  the 
solution  of  inhomogenous  boundary  value  problems  on  infinite  intervals. 


-10- 


*i^*0**WNK  'j- 


3.  The  Case  G(°>)  *  0  . 

Again  we  consider  Che  problem  (1.1),  (1.2),  (1.3)  but  we  drop  the  restriction 
G(“)  »  0  .  We  again  assume  that  B  is  a  r_  *  n-matr.'.- 
The  following  assumption  will  be  needed: 

( I )  The  problem 

(3.1)  y '  -  t°A( t)  y  -  0  1  <  t  <  « 

n  n 

(3.2)  Byh(1>  “  0 

(3.3)  yh  £  C< [!,«•] ) 

has  the  unique  solution  yh  »  0  .  This  guarantees  that  the  r  x  r  -  matrix 

(3.4)  BFIMI)  is  regular 

and  therefore  the  inhomogenoua  problem 

(3.5)  y'  -  t°A(t)y  -  taG(t)f(t),  1  <  t  <  » 

with  the  boundary  conditions  (3.2)  and  (3.3)  has  a  unique  solution  for  every  f  e  C([1,*>]). 
Moreover  we  restrict  the  eigenparameter  X  to  an  open  and  connected  set  flee  with 
0  e  ft  so  that  the  matrix  A<“>)  +  XG(“>)  for  left  has  no  eigenvalue  v(X)  on  the 
imaginary  axis  and  therefore  the  matrices  G+  and  G_  are  constant  for  X  e  ft  . 

de  Hoog  and  Weiss  (1980a)  proved  that  all  eigenvalues  X  of  (1.1),  (1.2),  (1.3)  which 
fulfill  X  e  ft  are  isolated  and  that  the  associated  spectral  subspaces  ate  finite 
dimensional.  Each  (generalized)  eigenfunction  y  associated  with  an  eigenvalue  X  e  ft 
fulfills  y(«)  «  0  .  The  spectrum  of  (1.1),  (1.2),  (1.3)  has  no  finite  limit  point  in 
ft  . 

Of  course  this  settles  the  case  G(“)  «  0  completely  because  then  ft  ■  C  holds  but 
the  compactness  arguments  in  Chapter  2  were  included  because  they  will  be  used  in  Chapter 
4  where  imaginary  eigenvalues  of  At*1)  will  be  admitted. 


-11- 


We  define  the  operator  V  slightly  different  to  chapter  2: 


(3.6) 


V  : 


Cfl((1,-1 >  ♦  CQ< [1,-J ) 


Vf 


where  CQ (  [  1 , 00 1  )  is  the  Banach  space  of  all  function  f  €  C([1,“))  which  fulfill 
f(°°)  *  0  and  y  is  the  solution  of  the  problem  (2.8).  (2.9),  (2.10).  Assumption  (I) 
makes  V  well-defined  on  CQ ( [ 1 )  and  (2.6),  (2.11)  guarantee  that  y(“)  -  0  if 
f<»)  -  0  . 

The  eigenvalue  problem  (1.1),  (1.2),  (1.3)  is  equivalent  to 

(3.7) 
with 

(3.8) 

because  all  generalized  eigenfunction  associated  with  X  e  (2  are  in  C  (  1 1  ,°°1  )  and 
because  X  =«  0  is  no  eigenvalue. 

Now  let  us  consider  a  fixed  eigenvalue  li  =  j  ,  X  £  I)  with  algebraic  multiplicity 
m  and  acent  B  .  The  spectral  projection  is  again  given  by 


Vf  -  pf  ,  f  e  C  ( [1,“]  ) 


M  =  .  X  e  ft 


(3.9) 


E  <=  E(y )  *  2mT  f  (2  -  V)-1dz  :  C0(  (1,-n  ♦  C0([1,“)  ) 


where  the  circle  T  centered  at  y  contains  no  other  eigenvalue  of  (3.7)  and  the  image  of 
T  under  the  mapping  X  »  ^  denoted  by  ^-(T)  is  in  ft  .  E  fulfills  (2.24),  (2.25). 

We  want  to  approximate  the  generalized  eigenpair  (X, Range  (E(y)))  by  a  sequence  of 
nearby  eigenpairs  of  (1.4),  (1.5),  (1.6). 

Therefore  we  define  the  operators  VT  for  T  sufficiently  large 

C( [1,T] )  *  C( (1,TJ ) 


(3.10) 


-  Vf 
T  T 


where  xT  solves  (2.36),  (2.37),  (2.38),  S(T)  is  independent  of  X  and  fulfills  (2.40), 


(2.41).  So  the  VT's  are  defined  properly  and  fulfill 


IV  f  1 . 4  <  const.  If_lr, 

T  T  [1 ,7}  T  [ 1 ,TJ 


(3.12)  'W'll.T)  "  C°n't,T  'V[1,T)  * 

So  each  vT  Is  compact  and  has  an  infinite  sequence  of  eigenvalues  y*  which  have  the 
only  accumulation  point  0  .  The  associated  spectral  suhspaces  are  finite  dimensional. 

It  is  therefore  clear  that  the  finite  interval  problems  (1.4),  (1.5),  (1.6)  can  not  be 
us.d  to  approximate  continuous  parts  of  the  spectrum  of  (1.1),  (1.2),  (1.3)  which  may  very 
well  exist  outside  of  i!  . 

He  define  the  restriction  operator 


C0<[1,-1  )  *  C([1,T)) 

f  *  “  *  I  f.  mi 

T  [1,T] 


Then  for  every  sequence  Tn  *  “  the  sequence  of  spaces  C([1,Tn))  form  a  discrete 
approximation  A(Cn(M,«J),  K  C([1,Tn)>,  rT  )  for  the  space  CQ([1,»])  in  the  sense  of 


Stummel  (1970). 


A  sequence  fT  e  C([1,T])  is  said  to  converge  to  an  element  f  e  CQ  ( [  1  ,«•]  )  , 
n 

denoted  by  fT  *  f  ,  if 
n 


I f„  -  r  fl 0  as  n 

T  T  [  1  ,T  ] 
n  n  n 


A  sequence  of  bounded  operators  in  C([1/TnJ)  ,  is  said  to  converge  to  a  bounded  operator 
on  C  ([!,•])  ,  again  denoted  by  Ay  ♦  A  ,  if  for  every  f  e  CQ  ( 1 1  ,«•]  >  and  for  every 


sequence  f 


fT  ♦  f  implies  fT  ♦  Af 


He  will  drop  the  subscript  n  mostly. 

Taking  a  fixed  z  *  0  in  the  resolvent  set  of  V  and  in  j  (ft)  we  investigate 
(z  -  vT)_1  .  Setting  uT  -  (z  -  vT)_1fT  for  an  arbitrary  fT  e  C((1,T])  we  easily  find 


that 


(3.16) 


UT  “  z  <yT  +  V 

where  YT  solves 

(3.17)  -  t0(A(t)  +  1  G(t))yT  -  -1  t“G(t)fT 

(3.18)  ByT(1)  -  0 

(3.19)  S(T)yT(T)  -  0  . 

Defining  F(-j)  as  the  matrix  which  transforms  A(“)  +  ^  G(“)  to  its  Jordan  canonical  form 
we  derive  from  de  Hoog  and  Weiss  (1980b)  that  (3.17),  (3.18),  (3.19)  is  uniquely  soluble 
for  T  sufficiently  large  if 

(3.20)  l(S(T)F(-l)G+)-,l  <  const  as  T  ♦  “ 
and  the  estimate 

(3.21)  l(z  -  V  )  ’l  ,  <  const  ( z) 

f ol lows  if  <3.20)  holds.  This  bound  is  uniform  in  z  e  ,  where  is  compact, 

0  l  Kj  and  -^(K^e  0  (see  Kreiss  (1972)). 

This  analysis  also  shows  that 

(3.22)  (z  -  VT)-1  ♦  (z  -  V)'1 
uniformly  for  z  e  1^  .  Therefore  (3.20)  guarantees  that 

(3.23)  inf  I P  -  UTf  *  0  for  t  ♦  - 

where  p  e  j(K)  are  the  eigenvalue  of  v  and  pT  are  the  eigenvalues  of 
VT  .  Moreover  the  spectral  projections  fulfill 

(3.24)  Et(U)  -  ~  /  (z  -  VT)~’dz  ♦  E(p) 

(3.25)  11m  ranMET(p))  >  rank(E(p)). 

The  sets  Range(ET(p))  form  a  discrete  approximation  A( Range(E(p ) ) ,  H  Range(ET (p ) ) ,rT) 

T 

for  Range(E(u)).  (see  Grigorieff  (1975)). 

In  order  to  make  sure  that  rank(ET(u))  “  rank(E(p))  for  T  sufficiently  large  it  is 
sufficient  to  show  that  the  sequence  ET ( u )  is  discretely  compact  (see  Stummel  (1971)) 
because  E(p)  has  finite  rank. 
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We  recall  that  the  sequence  of  bounded  operators  A^  In  C( [ 1 , Tn) )  Is  discretely 

n 

compact  if  for  every  bounded  sequence  fT  E  C( 1 1  *  T  ])  there  is  a  subsequence  f^  so 

n 

that  A^,  f  is  convergent  to  an  element  in  Cg([1,“])  . 

\  nk 

We  write 

(3.26)  ETrT  -  rTE  +  (ETrT  -  r^E)  iCQ  (  C  t  ,«I  )  -  C(f1,Tj)  . 
obviously 

(3.27)  ETrT  -  rTE  -  ^  /  <*  -  V’lV^  ~  V>(*  ”  V)'1d* 

holds. 


e‘  -  t  A(t)em  -  0  ,  1  <  t  <  T 

T  T 


For  an  arbitrary  f  e  CQ<  )  the  function  e?  -  (vTrT  ”  rTv>f  e  C([1,T])  is  the 

solution  of  the  problem 

(3.28) 

(3.29) 

(3.30)  SlTle^T)  -  -S(T)(Vf)(T)  . 

Proceeding  similarly  to  de  Hoog  and  Weiss  (1980b)  we  can  express  eT  explicitely. 


BeT(1)  -  0 


We  substitute  Fe  »  e_  where  F  is  as  in  (2.1)  and  get  the  problem 
T  T 


e^,  -  t°J(«)eT  +  t“(J(t)  -  J(-))eT 


(3.31 ) 

J(t),  J(“)  are  as  in  (2.11),  (2.2).  Now  we  write 

(3.32) 


W  +  V2  '  51  E  ^  '  «2  E  ^ 


where  e_  ,  e„  fulfill 
T  T 


(3.33) 


6  (  t) 
T 


j+ 

t _ "  ,  a+l  _a+1. 

exp(—  (t  -  T  ) 


+  (Ht(J  -  J(-))eT)(t) 


(3.34) 


e  (t) 
T 


+  (Ht(J  -  J(») ) e ” ) ( t ) 


where  H_  is  a  suitable  solution  operator  of  the  problem 


<3.35)  z'  »  t0J(t)z  +  t°g(t)  ,  1  <  t  <  T  ,  g  E  C((1,T])  . 

We  choose 

(3.36)  (HTg)(t)  -  (HgXt),  1  4  t  <  T 

with  H  defined  in  (2.4)  where  g  has  been  set  to 


(3.37) 


g(  t) 


g(t>  ,  1  <  t  t  T 

g(T)  ,  T  <  t  4  » 


Because  H  is  bounded  on  [5,"]  independently  of  6  we  get 

(3.38)  IH  (J  -  J(»))l,.  4  const.  »J  -  J(“)«r.  <  Vo 

T  [o,T]  [o,TJ  ‘ 

for  6,T  sufficiently  large.  The  operator 

(3.39)  I  -  (^(J  -  J<»))  :  C([S,T])  ♦  C([6,T)> 

is  invertible  and  e^,eT  e  C(f6,TJ)  are  uniquely  defined  and  can  be  continued  to  (1,T)  . 
Inserting  (3.32)  into  the  boundary  conditions  (3.29),  (3.30)  gives: 


T 

BFeT(1) 

BFeT(1) 

5i 

0 

“W 

T 

S(T)FeT(T) 

S(T )Fe  (T ) 
T  J 

-S(T)(Vf)(T) 

de  Hoog  and  Weiss  (1980a)  have  shown  that 

(3.41)  (a)  lim  e*(T>  -  G+  ,  (b)  He"  -  r^  1 . 1  T  ♦  0 

T  ' 

hold,  where  t  is  defined  in  (2.12). 

A  block  system  of  the  form 


(3.42) 


X  °° 


where  5,  C  are  quadratic  matrices  is  uniquely  soluble  if  and  only  if  B,(C  -  55  1  A  ) 

are  invertible  and  the  solution  is 
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-(C  -  DB  -1  A)”1  DB  -1 


B  ”'  +  B  _1  At C  -  DB  _1  A)  ’51  ' 


tc  -  Bi  ~1  a)"1 


-B  A(C  -  DB  A)'1 


The  off  diagonal  matrices  in  (3.40)  are  invertible,  their  inverses  are  bounded  as  T  ♦  *>, 
the  matrix  in  the  (2,2)  position  converges  to  0  as  T  *  “  and  the  matrix  in  the  (1,1) 
position  is  bounded,  and  therefore  the  system  is  invertible  for  T  sufficiently  large. 


(3.44)  lim  eT(1)  -  0 

holds  because  we  get  from  the  series  expansion  of  (3.33): 


co  a+i 

r  exp(— (W-T  )) 

,eT  *  l  0  J1 


<  CKJ-JI*) ) 


oo  a+i 

eXp(a+T(w-T  J) 


where  w(t)  -  t  has  been  set.  The  right  hand  side  of  this  inequality  can  be  estimated 


(3.44)(b)  o(k,T)  “  C  max  max  ( IJ(t)-J(<“)  »exp(-~T(ta+1-Ta  1 ) )  (Ta  ^t0  S*) 

i-1(1)k  tC(«,T] 

where  k  is  the  smallest  (positive)  real  part  of  the  eigenvalues  of  J*  and  k  is  the 
dimension  of  the  largest  Jordan  block  with  eigenvalue  with  real  part  <  .  Obviously 
<j(K,T)  *0  as  T  ♦  “  and  (3.44)  follows  by  continuation  to  [1,TJ  . 

We  get  from  (3.40),  (3.41),  (3.43),  (3.44) 


(3.45) (a) 


(3.45) (b) 


-  - ( ( S ( T ) FG+ ) ~ 1  +  o(T) )S(T) ( Vf ) (T ) 
5^  -  o(T)S(T)(Vf ) (T) 


holds.  For  g  e  CQ  (  [  1 ,«»]  )  we  therefore  get: 

(3.45) (eXV  r  -r  VHz-Vj"^  -  (-Fe+(^(T)FG  j"1  +  Fe”o|,T>  +  o(T) )S(T) ( V( z-V)-1 g) (T)  . 
Obviously  h  “  V(z  -  V)_1g  is  the  solution  to  the  problem 

(3.46)  h*  -  ta(A(t)  +  ^  G(t))h  -  ~  t°G(t)g(t) 

(3.47)  Bh( 1 )  -  0 

(3.48)  htC((|1,»)l 
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He  define 


(3.49)  Aft, A)  -  Aft)  +  AGIt)  for  t  611,“) 
and  the  family  of  operators  H ( X )  :  C([6,«>])  -*■  C( ;6,“)  ) 

(3.50)  (H(A) g)  (t)  -  <Mt,A)  /  P+  (A  )*_1  ( s ,  A )  sag(  s)  ds  +  $(t,A)  /  P_  (A  )$" '  (  s,  A )  s“g(  s)  ds , 

for  6  >  1  where 
(3.51  ) 


ta+1 

$(t,A)  -  exp(A(*,A)  — ) 


(3.52)  (a)  P+(A)  -  F(A)D+F_1(X),  (b)  P_(A)  -  F(A)D_F_1A) 

A  -  1 

hold.  Obviously  H(0)  -  FHF  for  H  as  in  (2.4).  The  projections  P  (A),  P  (A)  are 
holomorphic  in  compact  sets  K  c  0  (see  Kato  (1966)). 

Using  the  techniques  of  de  Hoog  and  Weiss  (1980a)  and  Markowich  (1980b)  we  conclude 


that 

(3.53)  IH(A)»[S  m]  <  Cf A ) 

where  C ( A )  is  independent  of  6  and  bounded  on  compact  set  K  <=  0  .  Moreover  it  is  an 
easy  exercise  to  show  that  (H( A ) g( • ,A ) ) ( t)  is  holomorphic  in  K  for  all  t  £  [6,»]  if 
g(t,A)  is  continuous  for  t  £  (1,°"I  and  holomorphic  for  all  t  in  A  e  K  . 

Proceeding  as  in  chapter  2  we  rewrite  (3.46) 

h'  =  t°A(“,-)h  +  t“(A(t,-)  -  A(»,~))h  +  \aG(t)g(t) 
z  z  z  z 


(3.54) 


B(t,~) 


and  get 

(3.55)  h ( t )  -  $(t. ■£)"_<■£);  +  ( H (~) B(  •  »~) h)  ( t )  +  ^(H(^)Gg)  ( t ) 

where  the  columns  of  the  holomorphic  n  *  r  matrix  W_(— )  span  the  range  of 
1  r- 

P  (— )  and  C  e  C  .  Because  of  (3.53)  there  is  a  fixed  6  so  that 

-  z 

(3.56)  IH(-j)B(*,'j)l[6  <  V2  for  all  z  £ 

where  S  is  the  closed  disk  contoured  by  f  . 


Setting 


( 3  •  57 )  ( a)  4_O.X)  -  (I  -  HU)B(«,X))'V*,X)»MA) 

(3.571(b)  4<Gf,X)  -  (i  -  H(X)B(.,X))_1H(X)Gg  . 

The  general  solution  of  (3.46),  (3.48)  on  (6,“]  is 

(3.58)  h(t)  -  i])  <t,-)C  +  -  i|i(Gg,-)(t) 

-  z  z  z 

\)i_,iKGg,X)  can  be  uniquely  extended  to  (1,®1  . 

Proceeding  as  in  (2.19)  using  (3. S3)  we  get 

(3.59) (a)  Ii)i_(T,-j)I  «  o(T)  uniformly  for  z  £  T 

and 

(3.591(b)  I5I  ®  0( Igl ^  )  uniformly  for  z  £  V 

can  be  concluded  as  in  Chapter  2  using  (3.56). 

(3.56)  and  the  uniform  convergence  of  the  series  expansion  of 

analyticity  of  (Gg ,— )  (T)  for  z  £  I  ■  Therefore 

z  )J 

<VTrT  ‘  rTV)(z-V)"1g  “  (-Fe^(S(T)PG+)‘1  +  Fe~  o(T)  + 

(3.60) 

+  O(T)  )  (-^S(T)t|((Gg,“)  (T) 
z  z 

—  1 

UT  =  (*  -  VT)  Fe^  is  the  solution  to  the  problems 

(3.611(a)  UT  ”  taA(t»~)u,r  “  ® 

(3.611(b)  Bu  ( 1 )  -  -BFe*(1) 

T  z  T 

(3.611(c)  S(T)u  (T)  -  ^S(T)Fe*(T)  . 

We  set  similarly  to  (3.32) 

(3.62  Ha)  uT(t)  -  e*(t,|){|  +  e‘(t,i)C* 

where 

(3.621(b)  e*  ( t , X )  >  4 ( t , X ) 4“ 1 (T, X ) w+( X )  + (H ( X ,T ) B( • , X ) e*< • , 

(3.621(c)  e”(t,X)  -  4  ( t ,  X  )W  _  ( X )  +  (H(  X  ,T )  B(  •  ,  X )  e^(  •  ,X ) )  ( 

hold.  The  columns  of  the  holomorphic  matrix  W+(X)  span  the  range 


(3.57)  assures  the 


+  o(T ) )  . 


X))(t) 

it) 

of  P+(X)  and 


H(A,T):  C([6,T])  ♦  C([6,T])  so  that 


(3.62)<d> 


H(A,T)g  J-  H(A  )g  ,  git) 


g(t)  , 


g(T>, 


6  <  t  <  T 


t  >  T 


holds.  Because  of  (3.53)  the  equations  (3.621(b),  (c)  are  uniquely  soluble  for  all 
z  e  if  6  is  sufficiently  large  and  the  analyticity  of  e*(t,^)  for  all  t  e  [1,T] 
follows  by  the  above  argument  and  by  continuation  from  [6,T]  to  [1,T]. 


Moreover  we  derive  as  in  (3.41)  from  de  Hoog  and  Weiss  (1980a)  (3.63) 

1 , 


(a)  lim  e*(T,-;)  -  W  (-) ,  (b)  le'(  1 


and  as  in  (3*44) 


T"'z'  ”+'z"  --  -T'-'I1  -  rT*-(,'I,’[1.T]  +  0 


(3.631(a)  lim  e*(1,-)  -  0 

_  T  2 

uniformly  for  2  e  r  . 

Inserting  into  the  boundary  condition  ( 3.61 ) (b) , ( c)  results  in  a  block  system  of  the 
form  (3.42)  and  using  (3.43)  gives 

( 3.64 )  ( a)  u  -  e*( •,-)<-  -(S(T)W  )_1S(T)FG  +  o(T) )  +  e‘(*,-r)°(T) 

i  1  Z  Z  +  Z  *  1  z 

uniformly  for  z  €  T  .  The  solvability  of  the  problem  (3.61)  follows  from  the 
invertibility  of  the  (analytic)  matrix  S(T)W+(^)  which  is  a  direct  consequence  of 
Assumption  (3.20).  Putting  all  together  gives 


(3.641(b) 


where 


( ( z  -  V  r’fV  r  -r  V)(z  -  v)_1g)  ( t)  = 

T  T  T  T 

=  e^(t,^)(S(T)W+(-j))_1S(T)i|i(Gg,^)(T)  + 

z 

+  e”(  (T,z)  +  (z-VT)_1Fe”(t)A2(T,z)  +  AjIT.z) 


I  A. ( t, z )  I  <  o(T) I gl  ,  for  i  =  1,2,3 

uniformly  for  z  e  T  . 

The  first  summand  on  the  right  hand  side  of  (3.641(b)  is  holomorphic  in  and 

therefore  its  contour  integral  along  T  vanishes. 


I 

i 


I 


1 
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Now  we  define  the  imbedding  operators 


(3.66)  (a)  i  iCUl.T))  ♦  C0<M.»]),  (b)  <iTfx><t> 


fT(t)  1  <t  <  T 


fT(T)  |  ,  t  >  T  . 


Obviously  rTiTfT  -  fT  and  [1#oo)  =  ,fTBC1,T]  holdB- 

Because  of  (3.64),  (3.65)  the  operator 

(3.67)  (z  -  VT)_1(VTrT  -  rTV)(z  -  V)-1iT  C( ( 1 ,T] )  ♦  C< [ 1 ,T] ) 

is  discretely  compact  : or  every  z  e  T  and  the  countour-  integral-operator  (3.27)  is  also 
discretely  compact  (see  Grigorieff  (1975).  Because  Range(E)  is  finite  dimensional 
rTEiT  is  discretely  compact  and  so  is  E^r^i^,  “  ET  and 

(3.68)  rank(ET<p))  »  rank(E(y)) 
for  T  sufficiently  large. 

Therefore  it  is  guaranteed  that  the  eigenvalue  u  •  ^  is  stable  with  regard  to  the 
VT‘s  (see  Grigorieff  (1975))  so  that  there  are  exactly  m  -  rank(E(uJ)  eigenvalues 
of  vT  which  converge  to  y  and  the  estimates 

(3.69)  gap(Range(ET),rT(Range(E)>)  <  const. «(VTrT  -  rxv  J  |  Range  ( E ) "  £  ,  T  ] 
and 


(3.70) 


max(  1 1 - 


max 

i 


|XT 


1 16) 


<  const 


1  ( Vr  - 
T  T 


rTV* I Range(E) *  \ 


1,TJ 


hold  (see  Grigorieff  (1975)). 

Under  the  assumption 

(a+1 )k_(X)+1 

(3.71)  A(  7  )  ,  G(  7  )  e  C  *  ([0,  J  )  ,  6  >  1 

where  k  (X)  is  largest  algebraic  multiplicity  of  an  eigenvalue  v(X)  of  the  matrix 
A(»)  +  XG(“)  with  a  negative  real  part,  the  estimates  (2.58),  (2.59),  where  v_  ■  \>_(X) 
is  now  the  largest  negative  real  part  of  the  eigenvalues  of  A(»)  +  XG(“)  ,  follow. 

However,  a  stronger  estimate  can  be  derived  by  proceeding  as  Osborn  (1975)  did  but 
without  carrying  out  the  last  estimates  which  lead  to  his  Theorems  1,2,3.  In  the  same  way 
the  estimates  given  by  Grigorieff  (1975)  can  be  changed.  We  get 


l  ' 
! 


1 
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{ 3 . 72 ) ( a)  gap! Range( Et  ) ,  rT  (Range! £)) )  <  const . I ( E^r^  -  rTE >  I RangetE ) 1 

(3 . 72  )  ( b)  max(|—  -  X|,  max  |X*  -  X|6)  <  const. HE  (V  r  -  r  V)  >  I  1 ' 

u  1-1  (1)m  T  T  TT  T  Range(E) 

T 

An  estimate  for  the  right  hand  side  of  (3.72) (a)  can  be  obtained  by  using  (3.641(b)  with 
g  e  Range (E)  and  (3.65): 

<v_(X)+£)  Q+1 

(3.731(a)  gap(  Range  (ET> ,  rT(Range(E) ) )  <  const.  o(T)exp( - -  T  ). 

In  a  similar  way  we  get 

1  3  (v_(X)+e) 

(3.731(b)  max(  |-s - X|,  maxlX^-  X|P)  <  const.  o(T)exp( - — -  T  )  . 

PT  1 

Therefore  the  standard  error  estimates  are  not  sharp  for  the  whole  class  of  problems. 
Retracing  the  history  of  the  o(T)  in  (3.73)  we  get 

(3.74)  o(T)  <  max(lA(T)  -  A(“)«,lG(T>-  G(»)l,  5( X  -  £,T)) 

where  o  is  defined  in  (3.44)(b)  and  £  is  small  when  the  radius  of  T  is  sufficiently  small. 

S (T)  can  be  chosen  independently  of  X  for  a  large  class  of  problems,  for  example  if 
G(— )  is  regular.  In  this  case  we  can  set  G(”)  -  I  because  this  always  can  be  achieved 
by  a  linear  transformation.  Then  F(X)  E  F  and  (3.20)  is  equal  to  (2.41).  ft  is  then  the 
strip  v  <  X  <  v+  where  v  is  the  largest  negative  and  v+  is  the  smallest  positive 
real  part  of  eigenvalues  of  A<“).  In  this  case  the  asymptotic  boundary  condition  (2.60) 
can  be  used. 

In  the  case  that  G(°“)  *  0  is  not  regular  S(T)  can  be  chosen  independently  of  X 
if  we  know  a  sufficiently  close  approximation  left  to  an  eigenvalue  X  of  (1.1),  (1.2), 
(1.3).  Then  we  rewrite  (1.1)  as 

(3.75)  y'  -  t“(A(t)  +  XG(t))y  -  7taG(t)y  ,  7  -  X  -  X  . 

A(t)  -  Ait)  +  XG( t) 

The  isolatedness  of  the  eigenvalues  X  £  ft  guarantees  that  the  problem 

(3.76)  (a)  y^  -  t°A(t)yh  -  0 

( 3.76 )  ( b)  Byh(1)  =  0 

( 3 . 76 )  ( c)  yh  e  C( [ 1 ,-]  > 
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Therefore  the 


has  only  the  trivial  solution  y.  =  0  if  A  is  sufficiently  close  to  X  . 

h 

above  theory  can  be  used  for  the  eigenvalue  problem  (3.75),  (1.2),  (1.3)  (with  7  as 
eigenparameter) .  The  change  consists  of  taking  X  instead  of  0  as  reference  point. 

We  set 

(3.77)  S  i  S(T)-  (G+)TF_1(X) 

where  the  superscript  T  denotes  transposition  and  F(X)  transforms  A(“)  +  AG(<“)  »  A(“) 
to  its  Jordan  canonical  form. 

Defining  F  «  F(A)  (2.40),  (2.41)  follows  inmediately. 

From  the  analysis  for  the  approximating  problems  (1.4),  (1.5),  (1.6)  it  follows  that 

it  is  sufficient  to  require  that  (3.20)  holds  locally  if  the  particular  eigenvalue  X  is 

to  be  calculated.  "Locally”  means  in  this  context  is  the  closed  set  oounded  by  the  contour 

— ( r )  defined  in  (3.9). 

U 

Since  the  family  of  projection  F(X)D+F  '(A)  is  holomorphic  in  K  c  ft  ,  K  compact 
there  is  a  nonsingular  r+  x  r+  matrix  T(A)  so  that 

(3.78)  w+(A)  «  F(X)G+T(X)  is  holomorphic  in  ft  . 

Therefore 

SF(X)G+T(X)  =  (G+)TF_1(X)(F(X)G+T(X)  +  0< |X  -  X| ))  -  T(X)  +  0( | A  -  X| ) 

holds  and 

(3.79)  (SF(X)G+)_1  -  T(  A  )T(  A)”1  +■  0(  I  X  -  X|)) 

for  X  sufficiently  close  to  X  ,  So  (3.20)  holds  locally  for  X  close  to  X  and  the 
asymptotic  boundary  condition  (3.77)  can  be  used  for  the  calculation  of  X  if  the  initial 
guess  X  is  sufficiently  close  to  X  . 

This  analysis  leads  to  the  idea  to  use  asymptotic  boundary  condition  which  depend  on 
the  eigenparameter  X  .  This  leads  even  in  the  case  that  the  'infinite'  problem  is  a 
linear  eigenvalue  problem,  to  nonlinear  approximating  'finite'  eigenvalue  problems  which, 
suggested  by  Keller  (1976)  have  been  successfully  used  in  computation  (see  Ng  and  Reid 
(1980)),  and  their  analysis  will  be  presented  in  a  subsequent  paper. 

However  for  many  important  fluid-dynamical  problems  it  is  possible  to  choose  simpler 
asymptotic  boundary  conditions.  An  example  is  presented  in  Chapter  5. 


4.  Imaginary  Eigenvalues  of  A(a°). 

We  are  now  going  to  neglect  the  crucial  restriction  that  all  eigenvalues  of  A(“) 
have  a  non-zero  real  part  but  we  will  require  a  sufficiently  fast  convergence  of  G(t) 
to  0  which  puts  us  back  into  the  compactness  argument  of  Chapter  2. 

We  assume  that 


(a+1)max(k  ,k  )+1 

(4.1)  A<  7  )  e  C  (  [0,  j  )  ),  6  >  1 

-<a+1)k  -e 

(4.2)  IG(t)l  <  const,  t  ,  e  >  0 

where  kQ  is  the  largest  algebraic  multiplicity  of  an  eigenvalue  of  A(“)  with  real  part 
zero  and  k_  is  defined  as  in  Chapter  2  . 

Markowich  ( 1980a, b)  has  shown  that  there  is  a  solution  operator  H  of  the 
inhcgnogenous  problem  (2.8)  which  fulfills 

(4.3)  l(Hf)(t)l  <  const,  t  E(int) If  I  j  1  ^  ,  1  <  j  <  n 

if  (4.2)  holds. 

Therefore,  if  the  homogenous  problems  (2.8),  (2.9),  (2.10)  has  only  the  trivial 
solution  y  =  0  then  the  operator  V  (see  (2.16))  is  well  defined  and  as  the  sum  of  a 
degenerate  and  a  compact  opeprator  it  is  compact  and  the  same  consideration  as  in  Chapter  2 
hold  for  the  eigenvalues  and  the  generalized  eigenvectors  except  the  decay  statements 
because  the  eigenvalues  with  real  part  zero  may  produce  solutions  which  are  asymptotically 
constant  or  which  decay  algebraically.  An  algorithm  which  determines  the  nature  of  the 
basic  solutions  under  the  assumption  (4.1)  is  given  in  Markowich  (1980b),  chapters  3  and  4. 

The  construction  of  the  supplementary  boundary  condition  S(T)xT(T)  «  0  for  the 
approximating  problems  (1.4),  (1.5),  (1.6)  now  relies  heavily  on  the  asymptotic  nature  of 
the  basic  solutions  and  is  explained  in  Markowich  (1980c),  Chapters  3  and  4.  The  matrix 
S(T)  constructed  in  the  mentioned  paper  takes  care  that  the  basic  solutions  which  are  in 
C([1,-J)  but  which  do  not  decay  sufficiently  fast,  are  dampened  by  the  multiplication 


with  S(T)  so  that  norm  convergence  of  the  operators  VT  defined  as  in  (2.35)  to  V 
results  (see  Markowich  (1980c),  Chapter  4). 


f 


Exponential  convergence  of  eigenvalues  and  spectral  subspaces  holds  if  all  (basic) 
solutions  of  the  problem: 

(4.3)  y'  -  taA(t)y  «  0 


(4.4) 


y  t  C{  (1,-) ) 


decay  exponentially 


5.  The  Orr-Sommerfeld  Equation . 


The  Orr-Sommerfeld  equation  (see  Ng  and  Reid  (1980))  governs  the  stability  of  laminar 

boundary  layers  in  the  parallel  flow  approximation 
2  2 

(5.1)  (2—.  -  a2  )2 *  -  [(U(z)  -  X)(iy  -  a2)*  -  U"(z)*]  -  0 

dz  dz 

a£R,a>0.  $(z)eia^x  is  the  disturbance  stream  function,  R  >  0  is  the 

Reynolds  number  and  U( z)  is  the  velocity  distribution  fulfilling 

(5.2)  U  e  C([0,®1)  ,  U(“)  »  1,  U"(»)  -  0  . 

The  boundary  conditions  for  the  Orr-Sommerfeld  problem  at  z  »  0  and  z  “  ~  are 

(5.3)  *10)  -  *'(0)  -  0 

(5.4)  *(»)  -  *'(«)  -  0  . 

This  problem  is  of  singular  perturbation  type  for  R  large,  but  we  disregard  that  and  just 
derive  appropriate  asymptotic  boundary  conditions. 

We  substitute 


and  get  the  problem 


0  10  0 
0  0  10 
0  0  0  1 


f,(z)  0  f2(z)  0 


1  y  =  X 


0  0  0  0 

0  0  0  0 

0  0  0  0 

a  0  b  0 


y  ,  0  <  z  <  » 


10  0  0 
0  10  0 


]y(0)  =  o 


y  e  c(  to,»] ) 


(a)  f^z)  »  -(a4  +  iaR(a2U( z)  +  U"(z))) 


(b)  f 2< *)  »  2a  +  iaRU(z) 


and 


(5.10)  (a)  a  «  ia2R  (b)  b  “  -ioR  . 

The  eigenvalues  of  A(°»)  are 

(5.11)  v  =  a  ,  v  -  (a2  +  iaR)  ^  ,  v  »  -a  ,  v  -  -(a2  +  ioR)^2 

12  3  4 

so  that  Re  v)(  Re  >0;  Re  ,  Re  v^  <0  and  the  eigenvalues  of  A(»)  +  AG  are 

(5.12)  VjU)  -  a,  v2(X)  -  (a2  +  iaR(l  -  A))1/2,  v^A)  »  -a,  v^A)  -  -(a2  +iaR(l  -  A))1/f2 

so  that  Re  v  (A),  Re  v„(A)  >  0  ;  Re  v  (A),  Re  v  (A)  <  0  for  all  A  e  C  holds. 

7  2  3  4 

Therefore  the  set  £1  defined  in  Chapter  3  is  the  whole  complex  plane.  We  get 


(5.131(a) 

and 

(5.131(b) 

so  that 

(5.14) 


J(“)  -  diag(v  ,v  ,v  ,v  )  for  A  *  1 
12  3  4 


for  A  *  1  ,  F(0)  ■  F 


(5.15) 


F(1)  - 


r 


1 

0 

1 

a 

1 

-a 

2 

a 

2a 

a‘ 

3 

a 

3a2 

-a 
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All  eigenvalues  X  are  Isolated  and  have  finite  algebraic  multiplicities 


holds  and  @  Is  the  ascent  of  X 


Computation  of  the  Orr-Sonmerfeld  problem  using  the  boundary  conditions  set  up  by 

(5.19)  can  be  found  In  Grosch  and  Orsag  (1977).  They  used  the  Blasius  velocity  profile 
2 

u(s)  "  1  +  0(e  ) ,  w  >  0  . 

-2aZ 

Their  numerical  experiments  indicate  that  the  order  of  convergence  is  e  in  the 

case  that  a  <  |Re  v^fX) |  <  1  and  X  has  (most  likely)  ascent  1  .  Checking  our  order 
formula  (3.73),  (3.79)  gives 


(5.22) 


o(Z)  <  max(e 


-wZ 


-wz  +(a-e)(r-Z) 


)  -  const. e 


(-a+e)z 


re[4,Z] 

and  the  order  of  convergence  the  theory  predicts  is  e 
spectral  subspaces. 


for  eigenvalues  and 
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